Possible interaction driven topological phases in (111) bilayers of LaNiC>3 
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We use the variational mean-field approach to systematically study the phase diagram of a bilayer 
heterostructure of the correlated transition metal oxide LaNiC>3, grown along the (111) direction. 
The Ni 3+ ions with d 7 (or ej) configuration form a buckled honeycomb lattice. We show that as 
a function of the strength of the on-site interactions, various topological phases emerge. In the 
presence of a reasonable size of the Hund's coupling, as the correlation is tuned from intermediate 
to strong, the following sequence of phases is found: (1) a Dirac half-semimetal phase, (2) a quantum 
anomalous Hall insulator (QAHI) phase with Chern number one, and (3) a ferromagnetic nematic 
phase breaking the lattice point group symmetry. The spin-orbit couplings and magnetism are both 
dynamically generated in the QAHI phase. 



PACS numbers: 

Introduction — Artificial transition metal oxide het- 
erostructures (TMOH) are becoming available owing to 
the recent development^ in the fields of oxide superlat- 
tices and oxide electronics. In particular, layered struc- 
tures of TMOH can now be prepared with atomic preci- 
sion, thus offering a high degree of control over important 
material properties, such as lattice constant, carrier con- 
centration, spin-orbit coupling, and correlation strength. 
The previous efforts on TMOH has been mainly focused 
on the (001) interface, where a rich variety of behavior 
emerges, such as superconductivity and magnetism(for 
a review, see Refi^). In addition, recent theoretical 
investigation^ pointed out that the bilayer TMOH grown 
along the (111) direction are promising materials realiz- 
ing various topological phases. 

The transition metal ions form a bulked honeycomb 
lattice in (111) bilayer structures (FigfTJ). Haldane first 
proposed that electrons hopping on a honeycomb lattice 
could realize the quantum Hall effect (QHE) in the ab- 
sence of Landau levels^, pointing out the possibility of 
nontrivial topology in simple band insulators. Such an 
insulator phase has been termed as quantum anomalous 
Hall insulator(QAHI). Its time-reversal symmetric gen- 
eralization, topological insulators, have attracted a lot 
of interest both theoretically and experimentally (for re- 
views see Refi^r— ). These insulators all feature a band 
gap driven by the spin-orbit coupling. In order to realize 
these physics experimentally, semimetallic materials are 
required. Honeycomb lattice is well-known to support 
semimetallic band structures, for instance in graphenc. 
Therefore the TMOH along the (111) direction is partic- 
ularly promising in searching for topological phases. 

The perovskite nickelates, RMO3 where R is a rare- 
earch atom, have demonstrated rich physics including 
metal-insulator transitions. One very interesting feature 
of these systems is a rather complex pattern of charge 
and spin orders (for a recent summary of experimental 




FIG. 1: (Color online) The structure of the (111) bilayer of 
RNi03. (a) The original perovskite structure, (b) The (111) 
bilayer of RNi03. Here, only Ni sites are shown, (c) Buckled 
honeycomb lattice formed in the (111) bilayer. 



progresses, see Ref J£). When R=La, the bulk compound 
remains metallic at all temperature. At low temperature 
it has a magnetic ordering pattern with an "up-up-down- 
down" spin configuration, coexisting with a "rock salt" 
type charge order. The charge order has been argued to 
be a by-product of the spin order based on symmetry 
considerations^. The magnetic ordering pattern can be 
visualized in the following wayii. When the cubic per- 
ovskite LaNiOa is viewed from the (111) direction, the Ni 
atoms form layers of triangular lattice. Each layer of the 
Ni atoms are ferromagnetically ordered. When these lay- 
ers are stacked along the (111) direction, the periodicity 
of the pattern is four: i.e., "up-up-down-down". Namely 
two adjacent layers are both spin up and the next two 
adjacent layers are spin down. After including the or- 
thorhombic distortion of the 3D lattice, this magnetic 
pattern is consistent with experimental observations^. 
Recently LaNiOa/LaAlOa superlatticcs along the (001) 
direction have been actively investigated experimentally 
and display unique quantum confinement effects^. 

Motivated by this interesting material, together with 
the recent experimental progresses of the growth of the 
(111) perovskite heterostructures (e.g. Refi^ 3 -), we study 
the possible quantum phases of LaNiOs bilayer TMOH 
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FIG. 2: (a) The band structure of the NN tight-binding model 
Eq.(HJ), each band is two- fold spin-degenerate. The lowest flat 
band is fully filled. This band structure can also be inter- 
preted as the one of the majority spin in the spin polarized 
DHSM phase found in our mean-field study, in which case 
unoccupied minority spin bands are not shown and each ma- 
jority band is non-degenerate. Fermi level is at the Dirac 
point, (b) The band structure of the spin polarized QAHI 
phase at Mo, y ,o = O.lt. Bands' Chern numbers are shown. 



grown on insulating substrates such as LaA103,LaSc03. 

Model — Because LaNiC>3 bulk material is metallic, 
a correlated itinerant electronic model would be a rea- 
sonable starting point. The bilayer forms a buckled 
honeycomb lattice (see Fig[lJ. The two e g orbitals, 
\d 3z 2_ r 2), \d x 2_ y 2) are not split in the trigonal environ- 
ment. Standard Slater-Koster constructional gives the fol- 
lowing nearest neighbor (NN) tight-binding model (t is 
the dda bond): 



H TB = - 2J z2 t fi» d ?, a ,<T d r t/ ,b,CT 
(r,r f ),cr ab 

t ( l -V3\ , t ( l V3 
1 o 
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Here r, r* label the positions of Ni, a, b label the orbital 
degrees of freedom, and a labels spin. 

The four bands (two sublattices and two orbitals) are 
shown in Fig 121 a), including two flat bands on the top 
and bottom of the spectra, quadratically touching with 
two dispersive bands in the middle at the T point. Sim- 
ilar to graphene, the two dispersive bands are linearly 
touching each other at the K and K' points, the corners 
of the Brillouin zone. For Ni 3+ ion with d 7 configuration, 
the t2 g is fully filled and the e g is 1/4 filled, so that the 
fermi level is positioned at the bottom quadratic band 
touching point in the non- interacting limit. 

It is well-known that the spin-orbit coupling in the e g 
orbitals is zero at the leading order due to the quenched 
angular momentum. The spin-orbit coupling can be in- 
troduced in the e g orbitals by higher order contributions 

in a trigonal environment. But because H$o = AL • S is 
weak for the 3d Ni 3+ ion: A ~ 80mcV— , simple estimate 
from the second order perturbation shows very small ef- 
fective spin-orbit coupling in the e g orbitals < lmeV. We 




FIG. 3: Results from the non-magnetic LDA+U cal- 
culations of the (a) LaA10 3 /LaNi0 3 /LaA10 3 and (b) 
LaSc0 3 /LaNi0 3 /LaSc0 3 TMOH with lattice relaxation of 
the LaNi0 3 bilayer. 



therefore do not include the atomic spin-orbit coupling 
in the tight-binding model. 

Since the 3d orbitals are quite confined in space, fur- 
ther neighbor hoppings are suppressed and this NN 
tight-binding model should be a rather faithful de- 
scription of the non-interacting electronic structure. 
In Fig|3] we present results from the first-principle 
GGA+U calculations of the LaA103/LaNi03/LaA103 
and LaSc03/LaNi03/LaSc03 bilayer TMOH. The cal- 
culations were performed by employing VASP code — 
in the context of density functional theory with the 
projector augmented wave (PAW) method 17 ' 18 for the 
atomic cores and the generalized gradient approxima- 
tion (GGA) 2& for exchange-correlation. The GGA+U 
method was used to treat the 3d electrons of Ni with the 
Hubbard on-site Coulomb interaction parameter U=7.0 
eV and J=0.65 eV— . With magnetism suppressed, these 
band structures are consistent with the NN tight-binding 
model, with t - 0.6eV for LaA10 3 /LaNi0 3 /LaA10 3 . 

The correlation on the Ni 3+ ions is expected to be in- 
termediate or strong. We consider the standard form of 
the on-site interactions: 

H I = Uy^ j niafUiai + J ^2 dl a , a dl bJ3 d la .pd ib . a 

i,a i.a<b 

+ U' n ^n lh + J Y (4a,t4a,l d *b,i d ^ + h.C.) 

z,o<6 i,a<b 

(2) 

U, U' arc intra-orbital and inter-orbital Coulomb repul- 
sions and J is the Hund's coupling. In our calculations 
below, for simplicity, we have set U' = U — 2 J, an equal- 
ity in rotational symmetric systems. 

The on-site U ~ 6 — 7eV has been used in LDA+U 
calculations for the nickelates (e.gj^l). However Ni oxides 
have strong charge-transfer effects^ to the oxygens. Our 
model H — Htb + Hi should be treated as an effective 
model. The values of U, J should be in a range such that 
the system has an intermediate to strong correlation. 

Symmetries — Here we only consider the inversion sym- 
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FIG. 4: (color online) The zero temperature mean-field phase 
diagram of H = H TB + Hi with U' = U - 2 J and J/U = 0.2. 
Vertical axis: the expectation values of the order parameters 
(defined in the figure, where di operator is given in the main 
text) in these phases. Numerical calculations were performed 
on a sample with I29 2 unit cells. 



metric case with the same substrate on both sides of 
the sample. The full lattice point group is D 3 d- Apart 
from the translational symmetry, the symmetry group of 
the system is D 3 d x SU(2) sp i n x TR where TR is time- 
reversal. The band touching points, both the quadratic 
ones and the linear ones, are protected by this group. 

Mean-field calculation — We have carried out the vari- 
ational mean- field study of the model H = Htb + Hi- 
After choosing J/U = 0.2, a reasonable ratio, the phase 
diagram as a function of U/t is systematically investi- 
gated. To be precise, we introduce the mean-field Hamil- 
tonian: 



H 



MF 



H 



TB 



(3) 



where i labels unit cells, and real numbers M a/ g 7 , with 
a, /3, 7 = 0, x, y, z, are the mean-field parameters, r, fi, a 
are all two-by-two Pauli matrices (the zeroth components 
are identity matrices) living in the sublatticc, orbital, 
and spin spaces respectively, di is the eight-component 
fermion operator including all these degrees of freedom. 
The mean- field ground state of Hmf, \MF), is used 
to minimize (MF\H\MF) numerically to determine the 
phase diagram. 

Because our model only has on-site interactions, only 
those mean-field parameters involving to. z are consid- 
ered. We have classified these order parameters accord- 
ing to the symmetry group D 3 d x SU{2) sp i n x TR, and 
find that each irreducible representation only shows up 
once. This means that the order parameters do not have 
to co-exist. 

The mean-field phase diagram for J/U = 0.2 at zero 
temperature is demonstrated in FigQ] As U/t is tuned 
from to 10, the phases and phase transitions are sum- 
marized in the following. When U — the chemical 



potential is at the quadratic band touching point at T 
point(see FigJ^a)). A spin nematic phase (SN) occurs at 
weak U/t. After a small region of coexisting with an un- 
saturated ferromagnctism, a first order transition drives 
the system into a fully polarized Dirac half-semimetal 
(DHSM) phase (see Figl2Ja)), where four bands of the 
the majority spin are half-filled and the fermi points are 
at the K,K' points. The spins remain fully polarized 
for larger U/t. Followed by a second order transition, a 
QAHI phase emerges(see FigJ5Jb) for its band structure). 
Finally after another first order transition, the system is 
in a fully polarized nematic phase (FPN). 

All spin orders are found to be colinear. For the discus- 
sion below, it is helpful to introduce a specific symmetry 
transformation in such colinear phases. Let us assume 
the order is along the S z direction. We define TR* to be 
a 180 degree spin-rotation (e.g. around the S x axis) send- 
ing S z — > —S z followed by a TR transformation. A usual 
colinear order respects TR* symmetry. In fact, a TR* 
symmetry-breaking immediately indicates a dynamically 
generated spin-orbit coupling. 

These phases are characterized by their order parame- 
ters and symmetry breaking. The SN phase (M 0jZiZ ^ 0) 
breaks D 3d , SU(2) spin and TR symmetries. The DHSM 
phase (M ,o,x ¥= 0) breaks both SU(2) spm and TR. In 
the QAHI phase (M ,o,, # 0,M 0)V ,o 0)&, SU(2) spin , 
TR and even TR* are broken. QAHI phase is the 
only phase in the phase diagram breaking TR* and 
has dynamically generated spin-orbit coupling. Finally 
the FPN phase (M 0fi , z + 0, M z , zfi + 0)24 breaks D 3d , 
SU(2) spin and TR. 

The QAHI phase, a band insulator with topologically 
protected chiral edge metallic modes, is also character- 
ized by a topological index — the Chern number or 
TKNN index^. The total Chern number of this phase is 
one (see Fig[^b)), which dictates quantized Hall conduc- 

tance o~ xy — — in the ground state. Such a dynamically 
generated QAHI in an SU(2) sp i„ symmetric Hamiltonian 
was proposed befor o 24 i 25 . 

We have found a dominant fully-polarized ferromag- 
netic order over the majority of the phase diagram. This 
tendency may be viewed as the residual of the bulk 
magnetic order "up-up-down-down" pattern, and can be 
qualitatively understood based on the large density of 
states from the flat band. Therefore we believe that it 
could be a reliable prediction of this mean-field investi- 
gation. The SN phase occurs at weak correlation and is 
unlikely to be realized. Note that the bandwidth of the 
bilayer system is substantially smaller than that of the 
bulk system due to coordination number reduction, and 
consequently the correlation in the bilayer should be even 
stronger than that of the bulk. This leads us to believe 
that only the DHSM, QAHI and FPN phases are within 
the reasonable regime of the real material. If DHSM or 
QAHI is found experimentally, it will be the first realiza- 
tion of such novel phases of matter. 

Concluding remarks — We have carried out a system- 
atic mean-field study of the phase diagram of LaNiC>3 
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bilayer TMOH, grown along the (111) direction. We 
hope that this study could encourage the experimental 
growths and characterizations of this system. Several in- 
teresting candidate quantum phases are found. Among 
them, the DHSM phase, similar to a spinless graphene, 
hosts symmetry-protected 2D Dirac cones. This phase 
has anomalous responses to an orbital magnetic field and 
can be detected by e.g. quantum oscillation experiments. 

Naively all the spin ordered phases could be destroyed 
at a finite temperature due to the Mermin- Wagner theo- 
rem. But the correlation length of the 0(3) nonlinear 
sigma model, a legitimate model describing the mag- 
netic order fluctuations, diverges exponentially at low 
temperature. This indicates that even an extremely 
small atomic spin-orbit coupling, which still exist in the 
LaNi03, could pin the order direction and support a 
rather high temperature phase. 

The tendency of developing a QAH gap in the DHSM 
phase can be qualitatively understood. At low energy, 
the only two possible gap terms in the DHSM are CDW 
(M Zi o,o) and QAH (M , y ,o)- CDW is disfavored by inter- 
orbital repulsion U' so the natural gapped phase contin- 
uously connected to the DHSM should be QAHI. 

The ferromagnetic order in the QAHI could form 
smooth textures — the skyrmions since ^(S 2 ) = Z. The 
dynamical generated spin-orbit coupling, which is an 
additional Z2 order parameter labeling the breaking 
of the TR* symmetry, also could form spatial domain 
walls. These topological objects could lead to novel 
physics. We point out that, similar to the quantum Hall 
ferromagnets2£, the skyrmion here is topologically bound 
with an electric charge and thus is a fcrmion. 

We have not discussed the possible Mott insulator 
phases with charge fluctuations completely suppressed. 
This possibility cannot be ignored particularly because 
the bandwidth of the bilayer is reduced significantly 



compared with the bulk LaNi03. Deep in the Mott 
regime, our model Hamiltonian is reduced to the Kugel- 
Khomskii-type2£ model whose leading terms favor ferro- 
magnetism on the mean-field level (e.g. see Eq.(2.7) in 
Re^&). If the full spin polarization persists in this regime, 
a simple t/U' expansion gives the NN model of or- 
bital fluctuations: H Mott = J£ ?£i (|U^ +iJ + Mp/^+jj + 

fip/xp,.), where A labels one sublattice, J = > 

and [f = (j, z ,[i a ' b = — \ll z T 2 This quantum 

model of pscudospins (i?, somewhat similar to the Ki- 
taev model^, has been used to describe multifcrroic lay- 
ered iron oxides^ and has not been solved in a controlled 
fashion. In a mcan-field-type study carried out here, ob- 
viously the anti-alignment of the orbital pseudospins is 
preferred which is exactly the characteristic of the FPN 
phase. However, quantum fluctuation could lead to ex- 
otic phases of matter and this forms a subject of future 
investigation. 

Finally we remark on some other possible phases at 
weak correlation. In fact, the quadratic band touch- 
ing point is known to be unstable towards interactions^. 
Apart from the nematic phase found here, this instabil- 
ity could lead to a quantum spin Hall insulator (where 
the spin-orbit coupling also comes from the spontaneous 
symmetry breaking), or a QAHI with Chern number two. 
Our mean-field calculation at J/U = has found these 
alternative phases at small couplings (U/t < 2). After we 
finished this manuscript, we have noted a related work^ 
focusing on these interesting phases. 
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